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We study quantum entanglement distribution on networks with lull-rank bi-partite mixed states 
linking qubits on nodes. In particular, we use entanglement swapping and purification to par- 
tially entangle widely separated nodes. The simplest method consists ol performing entanglement 
swappings along the shortest chain of links connecting the two nodes. However, we show that this 
method may be improved upon by choosing a protocol with a specific ordering of swappings and 
purifications. A priori, the design that produces optimal improvement is not clear. However, we 
parameterize the choices and find that the optimal values depend strongly on the desired measure 
of improvement. As an initial application, we apply the new improved protocols to the Erdos-Renyi 
network and obtain results including low density limits and an exact calculation of the average 
entanglement gained at the critical point. 

PACS numbers: 03.67.Bg, 64.60.ah 



I. INTRODUCTION 

The quantum repeater has been at the center of numer- 
ous studies addressing the distribution of entanglement 
over long distances, an essential prerequisite for many 
tasks in quantum information processing |TJ [2] . The cen- 
tral idea of the quantum repeater is to send one of a pair 
of entangled particles {e.g. a photon) across a series of 
links such that each link is short enough that the proba- 
bility of absorption is low and then to perform entangle- 
ment swappings at each node to further propagate the en- 
tanglement. However, the inevitable presence of noise in 
producing and transporting quantum states renders the 
straightforward application of repeaters hard in practice. 
This has resulted in the expenditure of a great deal of 
effort in designing distribution protocols. While the ini- 
tial repeater schemes involved a one-dimensional chain of 
nodes (containing qubits) connected by links, considering 
higher-dimensional networks of nodes and links has been 
a fruitful approach because the entanglement in neigh- 
boring links can be concentrated via purification. Given a 
fixed amount of entanglement per link that can be gener- 
ated in a particular laboratory setting, this concentration 
allows a useful amount of entanglement to be distributed 
over larger distances, which in turn allows protocols that 
consume entanglement to work over a larger distance. 
Previous work on entanglement distribution in dimen- 
sion greater than one considered either pure states [3HS] , 
or certain mixed states of rank two and three [6, 7 . Pure 
states were used in previous distribution studies because 
entanglement is better understood and easier to manipu- 
late in pure states than mixed states. However, realistic 
noise models imply that a bi-partite system of non-local 
components may only be prepared in a full-rank mixed 
state. In fact, it has been shown that long-range entan- 
glement on a cubic lattice of full-rank mixed states is 
possible [8]. 

On the other hand, while the bulk of work to date has 



been concerned with regular lattices, only a few works 
have treated random networks. The creation of entan- 
gled sub-networks on the Erdos-Renyi (ER) model [5] has 
been studied. Other workers investigated the effect of a 
particular transformation on several pure-state complex 
networks as well entanglement swapping with full-rank 
mixed states on the same networks [TOl E] . These stud- 
ies showed on one hand that transforming network hubs 
to rings via LOCC can enhance entanglement distribu- 
tion on a variety of pure-state networks and on the other 
hand how distribution is affected by the interplay be- 
tween the correlation length and the characteristic length 
of decay of fidelity under swapping. However, detailed 
studies of the application of of entanglement distribution 
protocols using both swapping and purification on com- 
plex networks have not yet been done. In the present 
work we address this deficit by introducing some natu- 
ral optimization problems in distributing entanglement 
along paths of full-rank mixed states. We find that the 
solutions to these problems yield surprisingly non-trivial 
results. Then, as a first application, we apply entangle- 
ment concentration protocols to the ER network. One of 
our main objectives is to understand when the use of the 
network connectivity offers an advantage for distributing 
entanglement between two nodes with respect to the sim- 
ple protocol in which entanglement is swapped along the 
shortest path connecting them. Note that this is always 
the case for a classical network: connectivity always helps 
in distributing classical information through a network. 
However, this may not be the case in the quantum regime, 
as quantum information cannot be cloned |12j . Indeed, 
we provide instances where the simplest direct protocol 
is better than the considered protocols using the network 
connectivity. 

In Sec. [Hj we give an overview of the models that we 
will examine. We shall consider the average concurrence 
on networks in which each link is initially a full-rank 
mixed state on two two-level systems, while a node is 
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FIG. 1: Part of a quantum network in its initial state. Small 
circles are qubits. Solid lines are bipartite states. Dashed 
circles enclose qubits within a node. Local operations may 
act on all qubits within a node. 



a 

(a) O 



1 ,'' 

-40 Q- 



O 



GO O- 



O 



-o 



o 



-o 
o 



(d) 



o 



-o 



FIG. 2: Entanglement swapping: (a) Before swapping a and 
7 are entangled and 8 and j3 are entangled, but systems cey 
and 8/3 are in a product state, (b) After swapping, systems 
a and /3 are entangled, while a/3 and 7<5 are in a product 
state. Purification: (c) Before purification, two entangled 
pairs (links), (d) after purification, one pair of nodes has 
higher entanglement than either of the original pairs. 



a local collection consisting of one party from each link 
terminating at that node. (See Fig. [I]) Unless stated ex- 
plicitly, when we speak of average concurrence we mean 
an average over both the outcomes of quantum measure- 
ments and the distribution of links for a given random 
network. We shall furthermore consider only two quan- 
tum operations for distributing entanglement (See Fig. 
21): 



• entanglement swapping, which probabilistically re- 
places a series of two links by a single link that 
bypasses the common node. The output link is in 
general less entangled than the input links. 

• purification, which essentially replaces two parallel 
links (i.e. sharing the same two nodes) by a single 
link that is more highly entangled than either input 
link. 

The main reason for the restriction to two operations is 
that many techniques that are successful on pure states, 
such as multi-partite techniques [13], are difficult, at best, 
to translate to networks of mixed states. However, these 
two operations naturally give rise to a rich set of protocols 
whose design is determined by the quantities that are to 
be optimized. 
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FIG. 3: Establishing entanglement between nodes A and B. 
(a) The shortest path Vab between A and B; the geometry of 
the path is irrelevant, so we represent it by a straight line with 
individual links not shown. Other paths connecting A and B 
are not shown, (b) The shortest path Vab (solid line with a 
dashed segment) between A and B. Between the endpoints 
of subpath 5 (dashed segment) there is an alternate path 
A (dotted line), (c) Similar to (b) with three subpaths and 
corresponding alternate paths. 



In Sections pI||III[ and |IV[ rather than designing a net- 
work for a particular task, we accept a given network of 
mixed states as a constraint. Our goal is then to identify 
and solve questions of design that arise in creating proto- 
cols to accomplish entanglement distribution. The main 
question is: when presented with the option either to swap 
or to purify, which is the better choice? For instance: 

• Single purification protocol. Consider the scenario 
in Fig.[3|D in which we want to entangle nodes A and 
B using the shortest connecting path Vab while 
making use of a neighboring path A. We proceed 
by swapping at all nodes on the two paths S and 
A to replace each of them by a single link, then 
purifying these two links, followed by performing 
swappings on all remaining nodes. In Sec. |IV C] we 
compute the ratio of path lengths ISI/I'PabI that 
produces the largest average entanglement between 
A and B, finding a value of approximately 0.37. 

• If the goal in the previous example is instead to 
achieve a positive probability of entangling A and 
B with minimal entanglement per link, then the 
optimal ratio of path lengths takes the value 1 + 
ln2/ln([v / 5 - l]/4) « 0.409. 

• Consider the scenario shown in Fig. [3j;, which we 
call the multiple purification protocol. Here, instead 
of purifying a single pair of subpaths, we purify 
n pairs. In Sec IV D we compute the minimum 



entanglement per link required to entangle A and 
B in the limit of large n for this protocol. 



In Sec. Ill we define specifically the direct and quantum 
strategics mentioned above. In Sec. |IV| we analyze the 
protocols in a more detailed and quantitative way. In 
Sec. |V| we apply the single purification protocol (SPP) 
mentioned above to a particular random network — the 
Erdos-Renyi (ER) network. We present results for short 
shortest paths and relatively impure states. We also com- 
pute the exact asymptotic concurrence of the SPP includ- 
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Werner state also has the advantage that it is created 
from any mixed state by removing the off-diagonal ele- 
ments via a depolarization process, a procedure that can 
be realized by local operations and classical communica- 
tion. 



FIG. 4: Quantum protocol. (a),(b), and (c) show the pro- 
gression of a quantum protocol: (a) Shortest path between A 
and B Vab with subpaths of lengths fei, fe, ni, n%. Dotted 
lines show alternate paths of lengths mi -I- n\ and m,2 + n% . 
In this paper, we shall always require that the collection of 
alternative paths and Vab be mutually disjoint, (b) After 
swapping and purifying subpaths. Each line segment now 
represents a single link with the labels giving the resulting 
Werner parameter, (c) After swapping all links. 



B. Concurrence as a measure of entanglement 

In this paper we use concurrence j!5j as a measure of 
useful entanglement in the system. For the Werner state 
the concurrence is given by 



C(x) = max{0, (3a;- l)/2}. 



(2) 



ing all shortest path lengths at the critical point of the 
model parameter. Finally, is Sec. |VI[ we address the ef- 
fects on the protocols of noise in the unitary operations 
and measurements. 



II. ELEMENTS OF THE MODEL 

We first introduce the networks that we shall consider. 
We then describe entanglement swapping and purifica- 
tion in more detail. 



A. Network and initial quantum states 

Consider the generic network of nodes and edges shown 
in Fig [l] With each edge of the network, we associate 
two two-level systems forming a bipartite system with 
states on C 4 . Thus, each node of degree k is occupied 
by k qubits. In the following, we shall consider states 
diagonal in the Bell basis 

| |4> Qb ) = ±= ( |0a) + (-1) 6 |la» : a, b e {0, 1} j . 

In particular, as the initial state on each edge, we choose 
the Werner state [T4] 

1 — x 

Pw {x) = x |$ooX$oo| + —7- U, (1) 

which has fidelity F = f ($ o| Pw{x) |$oo> = (3a; + l)/4. 
It can be shown that the Werner state is that it is entan- 
gled for x > 1/3 and separable otherwise. All protocols 
in this paper attempt to entangle two nodes by creating 
a Werner state on a pair of qubits, one from each node. 
Because of its simplicity, the Werner state serves as a 
standard form, allowing a clearer exposition of the dis- 
tribution protocols than does a Bell-diagonal state. The 



Because our task is to entangle any arbitrarily chosen 
pair of nodes, we define the average concurrence of the 
network 

C(x) = N{ ^_ 1) Y.^C{a,P), (3) 

where x is the parameter of the initial state, C(a,f3) = 
C{x a ^) is the concurrence of the state x a .p between a 
and ft after applying some protocol, and 7r Q ,/3 is the prob- 
ability that this protocol succeeds. Thus this definition 
depends on the choice of protocol and furthermore as- 
sumes that the concurrence of the resulting state between 
a and /3 is zero with probability 1 — ft a ,f}- Note that this 
average is over both pairs of nodes, as well as probability 
of the success of the protocol. 

In particular, we judge a particular protocol to be bet- 
ter than the direct protocol if it yields a higher concur- 
rence averaged over measurement outcomes. When ap- 
plied to our protocols on two-qubit Werner states, the 
concurrence has at least two advantages over other en- 
tanglement measures in this respect. Firstly, the concur- 
rence is the unique entanglement measure that is linear 
in x, which makes analysis easier. Secondly, the concur- 
rence provides the extremal comparison in the following 
sense. Most of the interesting entanglement measures 
are either convex (for instance, entanglement of forma- 
tion) or concave (for instance, logarithmic negativity). 
Suppose that a given protocol has higher average entan- 
glement than the direct protocol if concurrence is used 
as the measure. In Appendix [A] we show that the proto- 
col also has higher average entanglement, if any convex 
entanglement measure is used rather than concurrence. 
Conversely, if the protocol is worse than the direct when 
judged by concurrence, then it is also worse when judged 
by any other concave entanglement measure. Finally, we 
note that for the state the concurrence and the neg- 
ativity are identical. 
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C. Operations for distribution and concentration of 
entanglement 

1. Entanglement swapping 

In this section, we review entanglement swapping, and 
present the result of applying the operation to Werner 
states. Consider a state of four qubits a, /3, <5, 7, such that 
(a, 7) is an entangled pair and (5, /3) is an entangled pair, 
but systems 07 and 5/3 are in a product state. Entan- 
glement swapping is a sequence of quantum operations 
that transfers entanglement leaving (a, (3) entangled and 
(5,7) entangled. (See Fig. [2}) In the case of pure states, 
the optimal swapping is effected by measuring (7, 5) in 
the appropriate Bell basis, and then performing a cor- 
rective unitary on /3 depending on the outcome of this 
measurement |16j . with the result being either a maxi- 
mally, or a partially entangled state on In the 
latter case, swapping is usually understood to include an 
attempted singlet conversion on a/3, so that the result of 
the entire operation is to leave (a, /3) in either a maxi- 
mally entangled state (if successful) or a separable state 
(if unsuccessful.) Mathematically, we consider entangle- 
ment swapping to be a map from C 4 ®C 4 to C 4 , with the 
reduction of dimensions resulting from applying a partial 
trace over the system 7^. One can show that if the initial 
states on (a, 7) and (5, /3) are both maximally entangled 
(eg Bell states), then the swapping operation succeeds 
with probability one, assuming perfect operations. In the 
case of Bell-diagonal mixed states, we cannot know which 
of the four states we have drawn from the classical ensem- 
ble and so cannot unambiguously interpret the result of a 
measurement. The best we are able do is to assume that 
we have drawn the most probable state l^oo): and pro- 
ceed with swapping based on this assumption. However, 
if we have drawn a state other than |$oo)i we are un- 
fortunately increasing the average classical population of 
the remaining Bell states. More precisely, given two Bell- 
diagonal input states whose eigenvalues are (A, B,C, D) 
and (A' , B' ,C , D'), the un-normalized output state after 
swapping is 

(AA' + BB' + CC' + DD\ AB' + BA' + CD' + DC', 
AC" + BD' + CA' + DB', AD' + BC" + CB' + DA'). 

(4) 

Using ([I]) and Q it is easy to compute that performing 
swapping on two Werner states with parameters x and 
x' produces a Werner state with parameter xx' . That is, 



Pw(x) ® pw(x') i-> pw{xx'). 



2. Purification protocol 



(5) 



Purification protocols operate on a collection of bi- 
partite mixed states, producing a smaller number of bi- 
partite states of higher fidelity than the input states [T7] . 
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FIG. 5: x'(x,x) — x (solid curve), x — x\ ow , where a;i ow is 
determined by x'(x, x\ ov ) = x. (dashed curve). Both curves 
cross the x-axis at x = 1/3 and x = 1. 



We will use the Bennett-Brassard-Popescu-Schumacher- 
Smolin-Wootters (BBPSSW) purification protocol, intro- 
duced by Bennett et. al. |18j . which attempts to replace 
two input Werner states with parameters x\ and X2 by 
a single, more pure, Werner state — i.e. a state with pa- 
rameter x' satisfying x' > x\ and x' > X2- The parameter 
of the state resulting from this protocol is 



x'(xi,x 2 ) = 



Xl + X 2 + 4.T1X2 

3 + 3xiX2 



(6) 



with probability 



1 + xix 2 



while failure results in two separable (i.e. useless) states. 
Usually, in the study of purification protocols, one is con- 
cerned with the asymptotic limit of repeated purifica- 
tions. However, in the present case we are concerned with 
a single application of ^ . One common situation we en- 
counter below is purifying two states with x\ = x 2 = x. 
Another question is: Given a state x, what is the small- 
est value of a;i ow < x such that, when the states x and 
x\ ow are purified, the result is not worse than both of 
them; that is xi ow for which x'(x,x\ ow ) = x. In Fig. [5] 
we plot x'(x, x) — x and x — x\ ow v.s. x. We see that, 
in order for purification to be useful, X\ and X2 must not 
be too different, and also that purification is most useful 
for x « 0.7-0.8; Finally, we note that (|6) is increasing in 
both x\ and X2, a fact that we will use below. 

It is well known that the Deutsch-Ekert-Jozsa- 
Macchiavello-Popescu-Sanpera protocol (DEJMPS), in- 
troduced by Deutsch et. al. 1!) , yields states of higher 
fidelity than the BBPSSW protocol when performing re- 
peated purifications. In fact DEJMPS operates on two 
generic Bell-diagonal states, producing an output state 
that is also Bell-diagonal. When applied to two Werner 
states pw(xi), and pw(%2) both protocols yield states 
with the same fidelity and with the same probability of 
success. However, in general, only the output coefficient 
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of the | $00) component is the same for the two protocols 
with the remaining three coefficients differing between 
the protocols. In the case that the two input states are 
a Werner state and a general Bell-diagonal state with 
largest eigenvalue A, Q gives a state with concurrence 
[(4A- l)x- 1]/2. This result, together with the fact that 
the swapping Q is commutative and associative, imply 
that, for protocols using a single purification, the result- 
ing concurrence is the same whether we use BBPSSW or 
DEJMPS. On the other hand, direct calculation shows 
that applying Q to two states, each of which is the re- 
sult of purifying two Werner states, yields a state whose 
concurrence is improved with DEJMPS. In the present 
work, only the results in Sec. |IVD| are non-optimal in 
this sense. 

Finally, we mention that multiparticle recurrence |20j . 
and hashing |21j protocols have been shown to be more ef- 
ficient than protocols operating on two copies. Improve- 
ments have also been made by optimizing (in part by 
computer) over a large class of local unitaries [H] rather 
than using the unitaries employed in BBPSSW and DE- 
JMPS. These protocols may give better results, but they 
are more opaque conceptually and less amenable to anal- 
ysis. Furthermore, the gains shown in other contexts are 
rather modest. Thus, we do not consider these more 
complicated protocols here. 



III. ENTANGLEMENT DISTRIBUTION 
PROTOCOLS 

A. Direct strategy 

Our task is to entangle qubits on two selected nodes 
A, B of the network. The most naive approach is inspired 
directly by the quantum relay: perform repeated entan- 
glement swappings along the chain of links in the shortest 
path Vab connecting A and B using the procedure sum- 
marized in ([5]). (See Fig. [3^,.) For instance, swapping 
the first two links replaces p\$ — p2,3 — Pw{x) with one 
new link pi^ = pw(x 2 ). We then swap the resulting link 
pi j3 with p 3A = pw(x), yielding p M = p w (x 3 ), and so 
on. Thus, after swapping along the n interior nodes in 
Vab we obtain pab = Pw{x n ). We call this scheme the 
direct strategy. It is somehow analogous to a classical 
problem of transmission on a noisy network with trans- 
mission probability x on each link. However, the analogy 
is not perfect because the classical transmission proba- 
bility between A and B x n is in principle useful for any 
x > 0, whereas in our quantum network any Werner state 
with x < 1/3 is separable and thus useless as a resource 
for quantum information tasks. The concurrence of this 
direct strategy on a path of L links connecting nodes a 
and j3 is 



C Diroct (a,/3) d = C(x L ), 



B. Quantum strategies 

We may improve on the direct approach by using quan- 
tum mechanical operations to concentrate entanglement 
on the shortest path Vab connecting A and B. In partic- 
ular, we employ purification schemes to transfer entangle- 
ment from neighboring paths to subpaths of Vab- This 
follows the general idea of concentrating entanglement 
along a "backbone" that we used in previous work [3]. 
But, in the present setting, we must introduce new tech- 
niques because we are not trying to generate Bell pairs, 
and we must treat random neighborhoods of the back- 
bone. In what follows, these more complex strategies 
that exploit the network connectivity are called quan- 
tum, although it is clear that the direct protocol is also 
quantum. 



1. Swapping and purifying 

We begin by presenting the elementary combination of 
the purification and swapping protocols described above 
that we shall use in all of the protocols appearing below. 
Consider two paths, one of n links and the other of m + n 
links, with identical Werner states pw{x) on each link, 
for instance, paths S and A in Fig. [3j3. We first per- 
form entanglement swappings on each chain resulting in 
two states pw(x n ) and pw{x m+n ) which share the nodes 
at their endpoints. We then purify these two states to 
obtain a Werner state with parameter given by 



p x (n,m) 



x n _|_ gjn+n _|_ ^_ x 2n+r 



3 + 3x 2n+m 
the operation succeeding with probability 

, 1 + x 2n+m 
TT x (n,m) = . 



(7) 



(8) 



using C defined in (2). 



It is not difficult to prove that ^ only yields an im- 
provement over swapping alone (i.e. p x (n,m) > x n ) if 
m < n. 



2. Quantum strategy 

Here we present the class of protocols that we study 
in the remainder of the paper. In subsequent sections, 
we will study particular cases of this class of strategies. 
These strategies yield a higher average concurrence than 
the direct strategy. Referring to Fig. [3^, we say subpath 
for the segment S of Vab that we will purify. We say 
alternate path for a path A disjoint from Vab that we 
use to purify the subpath S. To entangle a pair of nodes 
A and B, the protocol is as follows. (See Fig. 

1. Identify the shortest path Vab between A and B 
of length L. Or, if there is more than one shortest 
path, choose one of them. 
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2. Identify a subpath Si of Tab with end nodes a\, 
b\ and length Tlx, such that there is an alternate 
path A\ of length mj 4- n\ with < mi < n-i 
joining ai, 61 that is edge-disjoint with "Pas- Note 
that we cannot have mi < 0, because this would 
imply, contrary to our assumption, that V ab is not 
a shortest path. 

3. Repeat step 2 zero or more times, as shown in 
Fig. [3J:, finding subpaths of lengths ni and m.; + n.i 
edge-disjoint from all previously identified paths. 
As depicted in Fig. [J|b), we now have a collection 
of subpaths of lengths and Wj + rii together with 
subpaths of V ab for which there is no sufficiently 
short alternate path. 

4. Perform entanglement swapping at each interior 
node on each of Si and Ai, effectively replacing 
each path of length I with a single Werner state 
Pw(x l ). 

5. Purify each pair of states that resulted from a swap- 
ping on each pair of paths (Si,Ai), This results in 
a new path connecting A and B as shown in Fig. 

Be). 

6. Swap along the new path connecting A and B to 
create a new Werner state between A and B with 
parameter 



L — n\ — ri2 - 



X 1 = p x (ni,mi)p x (n 2 , m 2 ) ■■■X 



In the following discussion, we find it useful to remove 
the length L from all quantities with the following change 
of variables. 



V 



L 



L 



(9) 



Note that y, a», 6, € [0, 1], and that a, and bi + aj are now 
the fractional lengths of the subpath and alternate path, 
respectively. The average concurrence of this quantum 
protocol is then written 

CQ M (a,/3) ^Yl^bAC ^-^ n^.Mj , 

(10) 

where i and j index the purifications in some arbitrary 
order, and we have used ^ and 

The choice of subpaths is not specified in the steps 
above, but is rather determined by the desired outcome. 
Below, we give explicit conditions on the choice of sub- 
paths for optimizing different quantities: maximum size 
of interval in initial fidelity for which QM protocol is 
better; n that gives minimum initial fidelity (minimum 
x) for which QM protocol gives positive concurrence; al- 
lowed values n near L for which QM protocol is better; 
n that yields the highest concurrence for fixed x. These 
protocols can, in principle, be applied to any network. 

The basic measure of success of the QM protocol given 
in steps 1-6 above is that it must give a better average 



concurrence than the direct approach. In the remainder 
of the paper we shall often be concerned with the increase 
in concurrence resulting from using a QM protocol. We 
denote this increase between nodes a and /3 connected 
by a shortest path of length L by 



AC(c*,/3) =C QM (a,(3) 
=C QM (a,/3) 



C micct (a,/3) 
C(y). 



(11) 



Likewise AC is AC(a,/?) averaged over a network. We 
call the interval in x for which the protocol is successful 
in this sense (plus some reasonable criteria) the good in- 
terval. The good interval is determined by the following 
criteria. 

• Each pair of subpath and alternative path must 
give an improvement in fidelity after purification. 
That is, 



Py(a,i,bi) > y a \ 



(12) 



if it succeeds. This requirement is necessary to 
avoid protocols which are advantageous, but would 
be even better if this particular purification were 
omitted. 

• For y < 1/3, the QM protocol must give a concur- 
rence greater than zero. That is, 



Y[p v (ai,bi) > -. 



(13) 



We call the root of the corresponding equality y* . 

• For y > 1/3, the average concurrence of the quan- 
tum protocol must be greater than the concurrence 
of the direct protocol. That is 



AC(a,/3) > 0. 



(14) 



We call the root of the corresponding equality y^. 
One can show that ( 12 ) and ( 13 ) give lower bounds 



on y, while (14) gives an upper bound on y. Phys- 
ically this can be seen as follows. If the quantum 
protocol gives positive concurrence for some value 
of y, then it will continue to do so for larger values 
of y (this also follows from the fact that (Jsj) is in- 
creasing in both arguments). At the upper bound, 
the effectiveness of the purification is decreasing 
with increasing y (as seen in Fig. [5| , but the prob- 
ability of success does not increase fast enough to 
make up for the decrease in the resulting Werner 
parameter. 

Finally, we note that the case in which the input pa- 
rameter y < 1/3 (that is x < (l/i) 1 ^) is especially inter- 
esting. In this case, the QM protocol is not only better on 
average, but is better in a stronger sense in that C T>1Tect 
vanishes for y < 1/3. 
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IV. ANALYSIS OF QM PROTOCOLS 

A. Generic form of constraints 

In this section, we present the constraints in a form 
that does not provide additional insight, but is useful for 
later calculations. The constraints (12), (13), and (14) 



determining the endpoints of the good interval are each 
of the form 



>o, 



(15) 




0.3 0.32 0.34 0.36 0.38 0.4 0.42 0.44 0.46 0.48 0.5 

y 



where {e^} is a relabeling of all of the a, and b' i7 and K 
and Cj i are some numbers depending on {di} the par- 
ticular constraint. We are not interested in the details 
of this formula, but we use it as a tool for calculating 
quantities appearing below. The end-points of the good 
interval are determined by the root, y* between and 1 

of f(y*,{di}) = 0. The end-point determined by (15) is The region inside the closed curves is where pp (y) "0 
thus given by and thus SPP is advantageous. The curve cutting through 

the closed curves is y a = 2y and maximizes ACf^ p (;y) with 
x * — y* T _ respect to a. 



FIG. 6: Roots of AC^ p b p (y) = from (20 1 as a function 
of y and a for various values of b. From the outermost to 
innermost curve the values of b are 0,0.01,0.07,0.11,0.135. 
Curves are determined from closed-form solutions a — a(y, b) 
with the roots for b = in particular given by (21 1 and (22 1. 



We denote by y the root of 

f(y,{di =0})=0. 



(16) 



The last expression is useful for computing perturbations 
around the solution of equations of constraint that are 
formulated such that all parameters vanish. 



B. Properties of constraints 

These properties hold for all protocols described by the 
six-step procedure above. 

• Consider for the moment just a single subpath of 
fractional length a and the alternate path of frac- 
tional length a + b. We ask which is better: swap- 
ping along the subpath and ignoring the alternate 
path, or swapping along each of them and purifying 
the result. The threshold at which purifying yields 
a Werner parameter equal to the input is described 
by Ap y (a, b) — p y (a,b) — y a = 0, which defines 
the threshold in each of the parameters a, 6, and 
y implicitly as a function of the other two. This 
equation is easiest to analyze if it is reparameter- 
ized as Ap y (a, ca) = 0, that is, by eliminating b via 
b = ca. The parameter c is also interesting because 
it gives the fractional excess length of the alternate 
path relative to the subpath. It is not difficult to 
prove that: 

i) The threshold c is given by c = c(y a ) where 



Hi) c(z) increases monotonically in z = y a , so that 
the threshold c increases(decreases) monotonically 
in y (a). 

iv) The difference in Werner parameter Ap y (a, ca) 
is maximized by c = for any fixed z, but is max- 
imized by non-trivial z for fixed c. For instance, 
Ap y (a, 0) maximized over z is approximately 0.05 



and is given by a root of 3z 4 + 8z 2 
numerical value z ss 0.69. 



8z + 1 with 



Because the map ([£]) from y* to x* is monotonic, 
the order of end-points of the good interval is pre- 
served as L varies. In fact, the intervals are com- 
pressed with increasing L. Thus, we only need to 
analyze the rescaled inequalities. 



c(z) = In (2[1 +4z- Zz 2 }- 1 ) hT 1 ^). 



(17) 



ii) c(z) takes valid values (i.e. non-negative and 
real) only on z G [1/3, 1] where we define c(l) by 
lim^i c(z) = 1. 



• The two roots determined by ( 13 ) and ( 14 ) coincide 
at y* = 1/3. This is because C7 Diroct (y) vanishes 
for y < 1/3, and increases continuously for y > 
1/3. Thus y — 1/3 is the threshold above which 
subtracting C lDlrect (y) from C^ M is necessary to 
evaluate whether the QM protocol is useful. 

• It can be proved that the largest absolute increase 
AC for the QM protocols compared to the di- 
rect, occurs at y = 1/3 for all protocols, that is, 

the largest y for which the direct protocol gives 
^Direct = q Thig ig ghown in g 



C. The single purification protocol 

We consider here the case of only a single purification 
(the single purification protocol [SPP] ) in which we iden- 
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0.32 0.34 0.36 0.38 



FIG. 7: Average concurrence as a function of scaled Werner 
parameter y. The solid curve is Co.409,0 which corresponds to 
the minimum (at a = 0.409) of the leftmost curve for b = 
in Fig. |6] The dashed curve is C Dircct . Note AC > for 
y* < y < The lower limit y* is determined by the largest 
value of y where C QM = 0. The upper limit y^ is where the 
curves coincide. 




FIG. 8: Concurrence from single purification QM protocol 
vs. a — n/L. For 6 = and y = 1/3. This corresponds to the 
dotted line in Fig. [5] To produce the maximum concurrence, 
the optimal fractional size of the su bpat h is a ~ 0.369. The 
curve is generated by setting b = in ( 20 1, and plotting C(y = 
1/3) vs. a. 



tify only a single subpath and alternative path pair. This 
situation is shown in Fig. [3]d. We analyze the protocol 
finding optimal values according to the most interesting 
metrics. In this case there is only one factor in each of 



the products in ( 13 ) and ( 14 1, while it is easy to see that 



(12 1 is redundant. Then (13 1 and (14) become 



,1+6 



Ay 



l+a+b 



,2a+b 



g(y) > o, 



where 



g{y) 



1 — y for y < 1 / 3 
5y — 1 for y > 1/3. 



(18) 



(19) 




FIG. 9: Region in ab-plane for which the single purification 
QM protocol is b ette r than direct (ACf^ p (t/) > 0), for several 
values of y. [See ( 20 1.] Curves are labeled with corresponding 



value of y. The good area lies below each curve. For y < 
1/3, curves with increasing y (dashed lines) entirely enclose 
previous areas. For y > 1/3 the same holds for decreasing y. 
The enclosed areas vanish at the smallest y* ~ 0.309 and the 
largest y£ — 1/2. The curves are determined from closed-form 
solutions b = b(y, a) of the roots of ( 20 1. 



and y > 1/3, y* and yj, respectively. Explicitly, the 
increase in concurrence gained from using SPP is 

AC a SPP (y) = \ {y b [V + y - (y a ~ 2y) 2 ] - g(y)} . 

(20) 

Inspecting ( 20 ), we see that for fixed y and independently 
of b, ACf p±J (y) is maximized for a solving (y a — 2y) 2 . If 
we further maximize over b and y, it is not hard to see 
that AC^ p b p (y) assumes a maximum value of 1/36 at 
b = 0, y '= 1/3 and a = (log(3) - log(2))/log(3). To 
further illustrate the behavior of |20| ), we consider the 
simplest case when a = n/L, b = 0, that is, the shortest 
path and the alternate path are of the same length. The 
roots of ( 20 ) solved for a are 



a(y) = 0, a(y) 
for y > 1/3, and 



log(4y - 1) 

log(y) 



(21) 



a(y) = 



log [2y ± 2 



log(y) 



(22) 



for y < 1/3. In particular, we see that the point where 
the roots p2| coincide gives the value of a representing 
the lowest lower bound y* on y and is given by a = 
1 + ln2/ln([v / 5- l]/4) « 0.409, with (See Fig. [i}) 



2/1 



= (Vs- l) /4« 0.309. 



(23) 



In accordance with the discussion above, we require 
y,a,b € (0,1). We call the roots of (|l8j) for y < 1/3 



Thus, this is the optimum value of a to allow the QM 
protocol to succeed with minimum initial fidelity. On the 
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other hand, the roots (21) coincide at the largest allowed 
value of y, 



(24) 



Inspecting (21 1 and (22) we also see that the largest good 



interval in y is obtained for purifying the shortest sub- 
path, i.e. as a — > 0. However, the improvement in con- 
currence also vanishes in this limit. (See Fig. [8]). Also 
note, as shown in Fig. |8j that the value of a that max- 
imizes the concurrence is different from the value that 
allows minimum initial fidelity as computed above. 

We now turn to the case b ^ (That is, alternate path 
is longer than subpath.) Swapping with a single purifi- 
cation is in every way worse than if b = 0. This follows 
from noting that the only effect on SPP of increasing b is 
to introduce a more weakly entangled state as one of the 
inputs to the purification. In particular, there is a value 
of a = n/L above which the QM scheme offers no im- 
provement for any value of y. One can further show that 
the maximum value of b allowing positive AC^ p b F (y) is 
b = log(7/6)/log(3) ps 0.14. The region in the a6-plane 
for which the single purification protocol yields an im- 
provement is shown in Fig [9} 



D. Multiple purifications 

Having analyzed the case in which we are allowed a 
single purification, we now turn our attention to the op- 
posite extreme of unlimited purifications. We partition a 
fraction a of the shortest path into n subpaths of equal 
length and purify each subpath with an alternate path 
of equal length. We ask how this protocol performs as 
n — > oo and find that the increase in concurrence tends 
to a limit, with a lower bound on y for which the protocol 
is good given by y = (l/3) 3 /' 3_a ^. 

Consider n subpaths Si of V ab of lengths a.; , not nec- 
essarily covering all of Tab, each of which has a cor- 
responding alternative path Ai also of length a, (See 
Fig I};.) We first swap along each subpath and alter- 
native path, then purify the resulting pairs. Finally, we 
swap along the all the remaining internal nodes. In this 
case we obtain from (13) the inequality for y*, 



y]l(l + 2y^)-l[(l + y^)>0. 



Likewise the inequality for y^ obtained from ( 14 1 is 



2"3yJ](l + 2^)-3 n n( 1 



-3y+l > 



In order to investigate the case of purifying many pairs of 
short paths, we choose the simplest case, setting dj = a 
for all i and a = a/n. That is, we consider purifica- 
tions on a fraction a of V ab , in which we purify n pairs 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

a 

FIG. 10: Constraints on y for multiple purifications with 
n = l,3,oo. The solid curves are y*. The dashed curves are 
y^ from the roots of (261. y* for n — > oo is taken from (25 1, 



while all other roots are determined numerically. The values 
of y between each pair of solid and dotted lines are the good 
interval, ie. those for which the QM scheme is better than 
direct for the corresponding value of n. 



of paths, with each path of rescaled length a/n. The 
inequality for y* is then 



(§)\(l + 2i,«) B -(l + y*y 



> 0. 



The limit of the solution of the corresponding equality as 

n — > oo is 



r = (i/3) 3 /( 3 ~ q ). 



(25) 



One can show that the inequality for y^ as n — > oo is 

3^ +1 -y Q -3y + l<0, (26) 

which we solve numerically. The results are presented in 
Fig. 10 together with the same curves for a single purifi- 
cation and three purifications. We saw that the value of 
a giving the minimum possible initial entanglement for 
a single purification is strictly between and 1. How- 
ever, for two or more purifications, the minimum is at 
a = 1 (this can easily be proven, as well). Thus, if the 
goal is for the protocol to work for the smallest possible 
initial entanglement, then performing purifications along 
the entire path is best in the present case. Also, both y* 
and y£ decrease with increasing n, with the lowest ini- 
tial entanglement possible giving non-zero concurrence 
y = (1/3) 3 / 2 w 0.19 for for a = 1 and n -> oo. This 
result demonstrates that for multiple purifications the 
best protocol performs as many purifications on short 
subpaths as possible, rather than fewer purifications on 
longer subpaths. In this sense, purifying before swap- 
ping is favorable. It seems very likely that the asymp- 
totic limit mentioned above is the best one can do (with 
our two allowed operations) without resorting to using 
previously purified links in further purifications. 
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E. Asymptotic form of constraints 

Here we consider the form of the generic inequalities 
of constraint ( 15 ) for large L, in order to find simple ex- 



pressions for the roots, which in turn give the endpoints 
of the interval where the quantum protocol is advanta- 
geous. We must take care, however, because we have 
some choices when taking this limit. We consider two 
different classes of limits. The first choice is one in which 
we ignore the rescaled equations so that L becomes large 
with rii fixed. In other words, we are holding the lengths 
of the subpaths constant as L becomes large. In this 
case, we find that the leading nontrivial term in the root 
is of second order in 1/L. The other choice is to let L 
become large with a t = rii/L constant. In this case, the 
rescaled equations are unchanged in the large L limit so 
we only have to look at the asymptotic form of the rescal- 
ing y = x L . Thus, for large L with held constant, the 
roots are given by x* s» 1 + In y* / L so that the inter- 
val between two constraints decays as 1/L. That is, the 
length of the good interval in x decreases as 1/L, 



Now we treat the case of holding m fixed. We proceed by 
first taking the small a% limit of the rescaled equations, 
followed by the large L limit of the inverse scaling x — 



y 



1/L 



An expansion of the LHS of (15) to first order in 



both y and a, gives 



y =y 



dyf(yAoj) 



where y is the root of ( [l6| . Replacing a, by n,*/Z and 
using (a + be) e = et £ + (b/a)e 2 + 0(e 3 ) we find to second 
order in 1/L 



x = y L = y 



L*yd y f(y,{0}) 

L2 J2, K J c j,oy c 



c 3,0 



(27) 



where ( 15 ) was used to compute the final line. Before 
proceeding to examples, we make two remarks on the 
expansions, i) for some values of the parameters rij, the 
numerator in (27) vanishes so the that the leading term 



in 1 /L in the length of the good interval is of order three. 
it) In some cases, we want to find the limit (27) for only 
a subset of {a^}, with the others held constant. In this 
case we simply remove some of the from the sums. 

There are several protocols in which these limiting 
cases are of interest. We mention two of them. Con- 
sider for example a network disordered in such a way that 
most shortest paths have nearly the same length (namely 
L). That is, if the shortest path between A and B is of 
length L, then the available subpaths and alternate paths 
are most probably of length near L. Furthermore, we re- 



writing |«S| = L — q and \A\ = L — r with L large and 
q and r fixed. Thus q and r represent small deviations 
in the length of available subpaths and alternate paths 
respectively. In this case, we define a — qj L, /3 = r/L , 
and make the substitutions a = 1 — a, b = —f3 in (|18j), 



and compare the result with (15) taking {di} — {a,/3} 
to find the parameters Cj j. Applying (27), we find the 
interval 



In 3 

1? 



2 ( 3 <L 



2r) < x < x 



In 3 

3L2 



(3g - 2r) 



Thus (assuming the roots are analytic in 1/L) only for 
q/r > 2/3, does a a good interval exist for large L. 

In the single and multiple purification schemes above, 
we saw that the optimal length of the alternate path A is 
the same as that of subpath S, that is m = 0. However, 
an alternate path of length exactly n will not be available 
in general. The lowest order fluctuation in the upper 
limit of the good interval, as m varies about 0, is studied 
by examining the small b limit, with the result 



y 



y- 



m\ny(—y + 4y 1+a — y 2a ) 
L 2 (y + 4(1 + a^+a - 2y 2a ) ' 



V. 



APPLICATION OF SINGLE PURIFICATION 
PROTOCOL TO ERDOS-RENYI MODEL 



consider the scenario in Sec. IV C of a single purification, 



We consider the Erdos-Renyi (ER) random graph [23l 
|2"1] because it is easier to analyze than more complicated 
random graphs and gives us insight into the behavior of 
the purification protocols on more complicated graphs. 
In particular, we want to compute the average concur- 
rence under the single purification protocol (SPP) of sec- 
tion [IVC] on the ER graph. The ER model is constructed 
as follows. Begin with the complete graph of N nodes 
and N(N — l)/2 edges and then delete each edge inde- 
pendently with probability 1 — p. Before proceeding, we 
simplify the notation below by introducing m! = n + m 
so that the alternate paths are of length m! . In the fol- 
lowing, we call ctl(p) the density of shortest paths of 
length L and r)L, n ,m'(p) the density of SPPs of the given 
parameters (that is, the fraction of pairs of nodes that 
admit this SPP). In general there is more than one pos- 
sible position for the subpath of length n along the SP of 
length L, and J7L inj m'(p) includes an average over these 
positions. The most important results in this section are 

• At low bond densities (small p) the density of 
SPPs characterized by L, n, m! is proportional to 
the product of the densities of shortest paths of 
length L and length m! . That is rjL.n,m' oc a^a m i. 
The constant of proportionality is determined by 
the number of positions for the subpath. 

• At high bond density ^L,n,m' ~ &l] that is, most 
subpaths have an available alternate path. 

• At the critical point Np — 1, and as N increases, 
all shortest paths are equally likely and the network 
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contains a number of each possible SPP of order 
1. As N becomes large and the Werner parameter 
is near 1, that is, 1 — x is small, the concurrence 
gained by applying all the SPPs is AC ~ AN~ 2 (l- 
x)~ 4 where A is a constant that is easily computed 
numerically. 

As we saw above, the SPP configurations can be par- 
tially characterized by the numbers {L,n,m!) giving the 
lengths of the shortest path, the subpath, and the alter- 
nate path, respectively. In order to compute the average 
concurrence for a particular value of Werner parameter 
x, we need to know the densities for various L, n, ml of 
the shortest paths admitting SPP that are beneficial for 
this value of x. It would greatly simplify understand- 
ing the protocol on complex networks if we could write 
the densities of the SPPs in terms of simpler and bet- 
ter known quantities, such as the distribution of shortest 
paths. To pursue the connection between these quanti- 
ties, we compute below the density of all SPPs on the ER 
network in the small p limit and see that in this limit, 
the density of an SPP characterized by L, n, ml is propor- 
tional to the product of the density of shortest paths of 
length L and the density of shortest paths of length ml ' . 
That is, r]L,n,m>{p) ~ g(L,n)a L {p)a m '(p). The factor g 
is discussed below. On the other hand, we argue that, as 
p — > 1, the density of SPPs with fixed L,n,m' is given 
simply by J7L,n,m'(p) ~ &l(jp)- (That is, nearly all short- 
est paths admit SPP). Between these two limits densities 
are more difficult to compute. One might expect similar 
behavior on other networks that have few connections 
(small p on the ER network), or many connections (large 
p), but we have not yet studied other networks in detail. 

A. Low Bond Density 



^3,2,2 




P 



FIG. 11: Density of single path purifications 773,2,2 for L = 
3, n = 2, ml = 2, on Erdos-Renyi graph with TV — 200. Points 
are MC data. Curve for smaller p is small the p expansion 
p 5 (TV — 2) (TV — 3) (TV — 4). Curve for larger p is the asymptotic 
formula (I — p 2 ) N ~ 2 (I — p) . The small and large p regions are 
shown in more detail using the same data in Figs. [T2"| and |13| 
Error bars are not visible on the scale of the plot. 
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In the limit of low bond density p, the numbers L, n, ml 
are enough to compute the density of the corresponding 
SPP. We take p to be small enough that two or more 
SPPs are unlikely to be available for a single pair of end 
point nodes A and B. The probability for the SPP con- 
figuration is 

(28) 

or T}L,n,m> (p) ~ g{L, n)p L+m ' N L+m ' ~ 2 for large TV, where 

,, . [ L — n + 1 for ml ^ n, .„„. 

9(L,n)={ ^ 29 

[L — n+ l)/2 for ml = n, 



is computed in Appendix |B} Similarly, we can show 

+ 0(p I+1 ), 



(N-L - 1)! 

or &l(p) ~ p L N L ~ 1 for large N. It follows that 
VL,n,m'(p) ~ g(Lin)a L {p)a m >(p) in this limit. 



FIG. 12: Density of single path purifications 773,2,2 for L = 
3, n = 2, ml = 2, on Erdos-Renyi graph with N = 200. Points 
are MC data (See Appendix ID]) Curve is the small p expan- 
sion p 5 (N - 2) (N - 3) (TV - 4)~ 



B. High Bond Density 

On the other hand, when p is large enough that a short- 
est path of length L is rare, then nearly all shortest paths 
are of length less than L. If the SP does not admit an 
SPP with subpath length n, then an edge-disjoint alter- 
nate path must be absent in all L — n+ 1 positions, which 
becomes rare with increasing p and L. It follows that 
nearly all shortest paths of length L will allow an SPP 
for all possible n and ml . Let us consider in particular 
L = 3, n = ml = 2. The density of such SPPs for small p 
is to lowest order in p 773,2,2 (p) = p 5 (TV-2)(iV-3)(TV-4), 
as shown in Figs. [TT] and [l2j For large p, the density 
of this SPP is nearly the density of shortest paths of 
length 3, which in turn is nearly the probability that 
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0.001 
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FIG. 13: Density of single path purifications 773,2,2 for L = 
3, n = 2, m! = 2, on Erdos-Renyi graph with N = 200. Points 
are MC data. Curve is the asymptotic formula (1— p 2 ) N ~ 2 (l — 
p). Error bars are not visible on the scale of the plot. 
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FIG. 15: Gain in concurrence AC over the simple method as 
a function of both p and x on the ER network. Monte Carlo 
results for TV = 200, considering all SPPs with L < 8. 
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FIG. 14: AC„/X m / /L (see (201 given by all possible single 
purifications with L < 7, vs. the initial Werner parameter x. 
Note that this figure is independent of the structure of the 
lattice. (Some purifications for L = 7 are omitted for clarity. 



the shortest path is not of length 1 or 2. It is easy to 
show (see Appendix [B]) that for all N and p, u\{p) = p 
and a 2 {p) = (1 - (1 - p 2 )^- 2 )^ - p). We then have 
?73,2,2 ~ \ — (J% — (J\ = {l-p 2 ) N ~ 2 {l-p) for large enough 
p, as shown in Figs. [TT] and [13] Thus, as we argued 
in the beginning of this section, to lowest order in p, 
~ cz,cr m /, but this no longer holds for large p 
where asymptotically T\L,n,m' 

« <j L . Finally, we address 
the consequences of these observations for the average 
concurrence. Figure [15] shows the average concurrence 
as defined in [3] on an ER network as a function of both 
input Werner parameter x and bond density p. This plot 
illustrates several features of the above analysis. The five 
concentrations of density correspond to L — 3,4,5,6,7 
(larger L were not computed in the MC calculations). 
For small p, longer SPs and SPPs are more prevalent, 
and these require larger x to be effective. On the other 
hand, for large enough p, most of the SPs are of length 1 



and 2, which do not admit SPPs. 



C. Concurrence at the critical point Np = 1 

In this section, we compute the asymptotic average im- 
provement in concurrence AC of the critical ER network 
for x near 1 and large N and find that 



AC 



A 



N 2 (l~x) 4 



(30) 



This expression is interpreted as follows. The factor of 
/V~ 2 is the probability that an SPP with any particular 
L,n,m and position of subpath will occur. The factors 
of 1/(1 — x) come from multiple SPPs contributing at 
one value of x: A fixed value of x gets contributions from 
SPPs with associated SPs of length L» 1/(1 — x); there 
are order L such shortest paths; order L different sub- 
paths (of length n) for each SP; order L alternate paths 
for each subpath; order L positions along the shortest 
path for the subpath. 



A few comments on the range of applicability of ( 30 ) 



are in order. In addition to requiring large N and small 
1 — x, we require that contributing paths not be too large 
so that the tree-like approximation remains valid. The 
most crude bound is that contributing paths be smaller 
than the radius (largest geodesic) of the network. At 
the critical point, there is a single cluster of size of order 
TV 2 / 3 with all next-largest clusters growing slower than 
any power. It has been proven recently |25j that the 
radius of the incipient giant cluster on the critical ER 
graph grows as N 1 ^ and furthermore (in distinction to 
the subcritical phase) the smaller clusters have smaller 
radii. Our numerical simulations show that the the ra- 
dius of the largest cluster is aiV 1 / 3 , with a approximately 
equal to 3. Using this radius as a bound on the valid 
range of L together with L ss 1/(1 — x) in (301, we find 
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that AC < 81 AN~ 2 / Z . Thus, wc sec that the advantage 
of single-path purification vanishes with increasing N at 
the critical point of the ER model. We expect similar 
behavior on other critical models as they will also have a 
broad distribution of very long paths. On the other hand, 
if we fix N — cp 2 , we get asymptotically a 2 = 1 — exp(— c) 
and 03 = exp(— c), in which case we expect SPP to con- 
tinue to show an advantage. 

The calculation of (30 1 proceeds as follows. It follows 



from (28) that 



,(p=l/N) = g(L, n)p 2 = g(L, n)/N 2 



(31) 



for large N. At this value of p the calculation of the av- 
erage concurrence is simplified in that the contributions 
from each path admitting SPP have the same dependence 
on N. Figure 14 shows the contributions to AC for in- 
dividual triples L,n,m', each of which is effective over 
a range of a;. At Np = 1, we are in the low density 
regime and only one SPP is likely to be present between 
any pair of vertices. Thus, for any value of x, the total 
contribution at Np = 1 is found by summing over the 
contributions for each triple L,n,m'. With increasing L, 
the density of SPPs with nearly the same proportions 
(that is, a and b) increases. Thus, although all these 
SPPs are equiprobable, as x increases the contributions 
come from increasingly large L with the number of over- 
lapping ranges increasing without limit as x — > 1. In fact, 
the average concurrence is 



AC(x) 



1 

N 2 



L 3 f(y) dL, 



(32) 



for large N and x near 1. Here y = x L and f(y), which 
accounts for the sum over n and to', is computed in Ap- 
pendix[Cj We define h(s) via h(—\n(y)) — f(y), and use 
— hix~l — x~e for x near 1. Then integrating (32) 
over L gives 



AC{x) 



1 

N 2 
1 



L 3 h{eL) dL 



5 /i(s) ds = 



N 2 e 4 



(33) 



The limits on the integral are determined by the lowest 
lower bound (23) and largest upper bound ( |24| on y (as 
shown in Fig!~[6). We have s* hi = - ln(yj!jand s z * = 
— ln(t/*) with A ~ 6.5 x 10~ 5 determined by numeric 
integration. In Fig. [17] we see that the asymptotic result 
(33) is approached rapidly with increasing L. 



VI. SPP WITH NOISY OPERATIONS 

Until now, we have considered only the ideal case of 
perfect operations and unlimited resources, with the only 
noise being that inherent in the Werner state. Even so, 
the SPP yields only small improvements in concurrence, 
with a maximum improvement of about 0.03. However, 




FIG. 16: h(s) appearing in (331. Inset is a semi-log plot of 



h(s) showing the function vanishing at S;* and as a power. 
The curve was determined by numeric integration. 



N 2 AC 




FIG. 17: Average concurrence scaled by iV v.s e = 1 — x 
on the ER network with Np = 1. The points are computed 
by summing all contributions of single purifications of short- 
est paths with L < 300. The solid line is the asymptotic 
result N 2 AC = Ae~ 4 . This plot contains no MC, but rather 
assumes iV is large enough that each SPP occurs with the 
probability given by the low density expansion (31 1. At the 



rightmost of the plot are the contributions from the shortest 
SPP with L = 3, n = 2, m = 0. 



because the maximum concurrence gain from a single suc- 
cessful purifcation using either BBPSSW or DEJMPS is 
only 0.05, purification schemes typically involve repeated 
purification of a large number of copies. Thus, the SPP is 
conceived as a first step to investigate possibilities of en- 
tanglement concentration on complex networks and per- 
haps as a building block in repeated purification schemes. 
However, noise can prevent even the small improvement 
in concurrence from a single purifcation so that that fur- 
ther purification is not possible. In this section we briefly 
consider the effect of imperfect unitaries and measure- 
ments on the SPP protocol. We employ a particular noise 
model for which the effects on the BBPSSW protocol and 
the swapping protocol were computed in reference pQ. 
Here, a noisy operation is modeled by a convex combina- 
tion of the perfect operation and a totally depolarizing 
channel that acts only on the same subspace as the per- 
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feet operation. An two-qubit operation on qubits 1 and 
2 with reliability pi is described by 

1 



Oi2P = P 2 Oi d 2 Ca V + 



P2 



tri 2 {p} ® I12, 



(34) 



with a similar definition for a single-qubit operator of 
reliability p\. An imperfect operation on a single qubit 
with reliability p\ is described in an analagous way. The 
imperfect measurement in the computational basis of a 
single-qubit is described by the POVM 



PS 



77|OXO| 

v\m\ 



(i 
(i 



•»?)|iXi| 

■ »y) |0X0[ , 



(35) 
(36) 



which is a projective measurement only when the param- 
eter 77 is unity. In pQ, the effects of theses noisy opera- 
tions on each the BBPSSW or DEJMPS and swapping 
protocols was computed. Although the DEJMPS proto- 
col was reported to be much more robust against noise, 
it is also much less amenable to analysis and was thus 
treated numerically. Here, we present simple closed-form 
results using the BBPSSW protocol. Purification of two 
input states pw{x) yield a state pw{x') where 



(2x + 4x 2 )(l - 5) 
3(l + a)+3x 2 (l-2<5)' 



(37) 



with <5 = 277(1 — 77) and a = (1 — p^)/p\- The probability 
of success is {l+a+x 2 (l — 2S)]/2. It is evident that, when 
a = 5 = 0, these reduce to the result for perfect opera- 
tions given by ([6]). We assume that swapping with noisy 
operations a chain of n Werner states each of parameter 
x produces a state 



(38) 



with parameter c > 1, which allows us to track sepa- 
rately the effects of noise from swapping and purification. 
If we further assume that the error model for unitaries 
and measurements described above applies to swapping 
as well, then reference [T| gives 

c = 7T-9 TT- (39) 

PlP 2 (477 2 - 1) 

The rescaled Werner parameter for a geodesic of length 
L is now 



y 



Note that y must now satisfy y < (l/c) L rather than 
y < 1 as in the case of perfect operations. For simplicity, 
we restrict our attention to the case b = in which the 
alternate path has optimal length. The average concur- 
rence (11) with noisy swapping and purification is then 



1 j 4c 2 (l-<5) 




-a- c 2 (l - 25) 



V 2 - g(cy) - 2cSy 



2(1-5) 
1-25 "' 



(40) 
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FIG. 18: Maximum value of AC with noisy operations, as a 
function of P2 and 77, normalized to 1 for perfect operations. 



where 



. I 1 — w for 77j < 1/3 
q(w) = < 

yy ' [4w-l for w > 1/3. 



(41) 



To simplify the analysis further, we note that for fixed y, 
ACf pp (y) obtains its maximum value at a = a max (y) for 



which the squared expression containing y a in ( 40 ) van- 



ishes. Furthermore, the maximum over y is obtained for 
V = 2/max = V(3c). This maximum average concurrence 



A C ^ = iJM-l(l + 2i) 



4 I 9(1 - 25) 



(42) 



which is independent of c, is plotted in Fig. |18| 7j.g. p 2 and 
77, normalized to the value for perfect operations 1/36. 



We see that at y 



!Jn 



and a = 



c (2/max) for er- 



rors of a couple percent, SSP yields improvements of the 
same order as for perfect operations. This is consistent 
with the sensitivity of the BBPSSW protocol to noise. 
In fact, for y = y max and a = a max (y max ) the values of 
the noise parameters for which AC^ pp (y) vanishes are 
exactly those for which x = x' in p7]). 



VII. CONCLUSION 

We have introduced and solved optimization problems 
resulting from the interplay between entanglement distri- 
bution and concentration. Already for simple protocols, 
the optimal choice of parameters is non-trivial and de- 
pends strongly on the quantity to be optimized. 

There are many unexplored questions still to be ad- 
dressed. For instance, our approach deals with a static 
initial network and searches for a protocol with no con- 
sideration of dynamics. 
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Appendix A: Entanglement on two-qubit Werner 
States 



Appendix B: Shortest paths and SPPs on the ER 
network 



In this appendix we show that the concurrence is 
the extremal entanglement measure when comparing the 
quantum to direct protocols. Below, we show that en- 
tanglement measures can be parameterize by C, i.e. 
E = E{C). We label two values E t = L^C,) with 
i — 1,2. The condition (10) has the form 

(Al) 



In this appendix we compute the density of shortest 
paths and paths admitting single purification protocol. 



pC 2 >d, 



with < p < 1. We will show that ( Al I implies pE 2 > E\ 



if and only if E(C)/C is non-decreasing, a condition sat- 
isfied by all convex entanglement measures. Thus, if the 
quantum protocol is advantageous according to concur- 
rence, then it is advantageous according to all convex 
entanglement measures. 

In this paper, an entanglement measure is a function E 
from density operators to [0, 1] that satisfies the following 
conditions 

i) If p is separable then E(p) = 0. 

ii) E(Be\\ state) = 1. 

Entanglement measures usually satisfy one or both of 
two other properties that will not concern us here: LOCC 
cannot increase the expectation value of the entangle- 
ment, and for pure states E reduces to the entropy of 
entanglement |26j . Many useful entanglement measures 
are convex, that is, 
Hi) for positive pi and '^2 i Pi = 1, 



PiE(pi) >E \J2p*P 



Condition i implies that entanglement measures must 
vanish for x < 1/3, so we need not concern ourselves 
with these states. The concurrence (I2J) is an invertible 
linear function for states with x > 1/3, so they can be 
parameterize by C rather than x, with C € [0,1], and 
we write E — E(C). Because the eigenvalues of the 
Werner state ([!]) are linear in C, the set of states with 
x > 1/3 is closed under convex combinations, so that Hi 
implies E(C) is convex. Similar statements can be made 
about concave functions. We now show that (Al) im- 



plies pE 2 > E 2 if and only if E{C)/C is non-decreasing. 
Clearly, \Al\ is equivalent to p £ [C\ /C 2 ,l]. In the 
worst case, we must then have (Ci/C 2 )E 2 > Ex, that 
is, E(C)/C is non-decreasing. That the converse is true 
can be shown with similar arguments. Furthermore, it is 
easy to show that for all convex entanglement measures 
E(C)/C is non-decreasing. A similar argument shows 
that pE 2 < Ei implies pC 2 < C\ if and only if E(C)/C 
is non-increasing, a condition that is satisfied by all con- 
cave entanglement measures. Finally, it is worth not- 
ing that we can make sharper statements. For instance, 
E a {C) = (C + AC 2 - C 4 )/4 is neither convex nor con- 
cave, yet E a (C)/C is increasing. It can be shown that 
the inverse of E a , E™ V (C), is also an entanglement mea- 
sure that is neither convex nor concave, but E ln ~ v (C)/C 
is decreasing. 



1. shortest paths 



The density of shortest paths of length 1 is obviously 
a i{p) — P- To compute cr 2 (p), consider the possible 
path of length 2 between vertices v a and Vb that passes 
through v c . This path is absent with probability 1 — p 2 . 
Note that the collection of paths for each of the N — 2 
possible v c together with the possible path of length 1 
are mutually edge-disjoint. Thus the probability that 
there is a path of length 2, but none of length one is 
a 2 (p) = (1 - (1 - p 2 ) N - 2 )(l - p). We did not com- 
pute er 3 (p), which would be more difficult because in- 
dependence is no longer present. However, we can say 
something about the case of large and small p. For 
large p, a(p) l+i /< j (p)l vanishes with increasing p so that 
ctl(p) ~ 1 — cl-i(p) — vl- 2 {p) ■ ■ • > which allows us to 
compute the asymptotic form of 03 (p). On the other 
hand, for small p, the probability of more than one path 
Between v a and vi, of length L becomes negligible and <jl 
is then the sum of the probabilities for each possible SP. 
The number of ordered choices of intermediate vertices 
for the SP is (A-2)(A-3) . . . (N-L- 1), each of which 
corresponds to a path with at least one unique edge. So, 
for small p, <r L {p) = p L (N-2){N-3) . . . (N-L-l). Our 
Monte Carlo calculations show that this approximation 
holds for the case Np = 1. 



2. SPPs 



The computation of density of SPPs is similar to that 
of SPs. The number of edges present in an SPP with 
shortest path of length L, subpath of length n, and al- 
ternate path of length m! is L + m' . Likewise the number 
of intermediate vertices is L+m' — 2. However, only in the 
case n = m' , the possible permutations of the L + m' — 2 
vertices can be partitioned into pairs, in which each mem- 
ber of the pair defines an SPP including exactly the same 
edges. In each pair, the interior vertices of the subpath 
and the alternate path are swapped, while the remaining 
vertices are unchanged. This is the origin of the factor of 
1/2 in ( 29 ). Finally, the number of ways that the subpath 
can be placed along a path of length L is L — n + 1, from 
which we arrive at ( 29 1 . Monte Carlo data supporting 



this expression is shown in Fig. 19 
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FIG. 19: r/L,n,n/g(L,n) vs. N at the critical point Np = 1 
for eleven pairs of L and n , fro m 3,2, through 7,3. The solid 
line is TV -2 as predicted by ( 31 1 . Points are MC data obtained 



by generating ER network samples and counting the number 
of SPPs. 



Appendix C: Density of contributions to 
concurrence at the critical point on the ER network 



Here we compute f(y) appearing in (32). In the 
following we let p = 1/N. We write ACi )Tljm (x) = 
max[C^ m (a;) — C^^^x), 0] for the average increase 
in concurrence obtained from purifying an SPP with pa- 
rameters L, n, m between two vertices. (The average is 
over quantum outcomes, distribution of ER networks, 
and possible positions of the subpath.) Here we use m 
and m! = n + m. As above, we have b = m/L which takes 
values between and 1. Then the contribution to the av- 
erage increase in concurrence between a pair of vertices 
for a fixed value of L is 

L-l n-l 



AC L X! VL,n.,ra'{p)AC L ,n,m{x) 

^ L-l n-l 

= jp Y Y 9(L,n)AC L ^ m (x) 



n—2 ??i—0 

n=0 
n-l 

+ AC L<n>m (x) 

777—1 

l-l/L 

- w y (i-o+Vi) 

a=2/L 
a-l/L 

b=l/L 



-AC Ltrl fl(x) 



where we used (29), and a — n/L, b — m/L take discrete 
values. This expression holds for large N, but for all 
x and L. Now for large L we replace the sum with an 
integral and find 



ACr 



L 



2 rl 



da (1 — a) 



1 



-AC*i, aj0 (x) 



L f dbAC ha , b (x) 
Jo 



L :i 



L'> 



v2 i (l-a)AC ha Ay)dadb=—f(y). 



In the last line, we discarded the term that is of order L 2 
and kept the term of order L 3 . The final expression is 
valid for all x, but for any fixed value of x, ACl vanishes 
with increasing L because y = x L . In practice, we com- 
pute f(y) numerically by integrating ( 10 1 over values a 
and b for which AC is positive as shown in Fig [9] The 
integrand could be done analytically, but the boundaries 
in Fig [9] are determined via numerical roots in any case. 



Appendix D: Monte Carlo computations 



We used a modified version of the C language library 
igraph 27 for Monte Carlo calculations. In particular, 
we replaced the calls to the system random number gen- 
erator with the Mersenne twister [55] generator. The 
number of trials for computing statistics varies greatly 
with model parameters from a few tens to 10 9 . Error 
bars for quantities such as the number of SPs and SPPs 
were computed from the sample variance in the mean 
number of events per network as a/^/Wt where n t is the 
number of trials. In a few instances, for the largest er- 
ror bars, n t is not large and the error bars are thus not 
accurate. 
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